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Abstract 



We evaluate, within one-loop chiral perturbation theory, the two-pion ex- 
change diagrams with single and double delta-isobar excitation contributing 
to elastic NN-scattering. We find that virtual A-excitation processes (in the 
static limit) produce the correct amount of isoscalar central attraction as 
needed in the peripheral partial waves with L > 3. Furthermore we compute 
the two-loop diagrams involving the 7T7r-interaction (so-called correlated 2ir- 
exchange). Contrary to common believe these processes lead to negligibly 
small and repulsive corrections to the NN-potential. The exchange of vector 
mesons (p, uj) turns out to be important for the F-wave phase shifts above 
Tiab = 150 MeV. Without adjustable parameters we are able to reproduce 
the empirical NN phase shifts up to 350 MeV for L > 3 and up to about 
(50-80) MeV for the D-waves. This is therefore the characteristic window 
in which the NN-interaction is basically governed by chiral symmetry. Not 
surprisingly, the lower partial waves require non-perturbative methods and 
additional short-distance parametrizations of the NN-dynamics. 



*Work supported in part by BMBF. 



1 



I. INTRODUCTION AND SUMMARY 



In a recent paper Jl| we have developed a framework for the application of chiral 
perturbation theory to low energy elastic nucleon-nucleon scattering. In that framework 
a systematic expansion of the NN T-matrix in powers of small external momenta and 
quark masses is performed by evaluating tree and loop diagrams with vertices taken from 
an effective chiral Lagrangian. The latter is an efficient tool to implement the (chiral) 
symmetry constraints on the dynamics of pions (the Goldstone bosons of spontaneous 
chiral symmetry breaking in QCD). Since such an approach to NN-scattering focuses on 
the dynamics of pions it is expected to work within the kinematical domain where the NN 
force is dominated by one- and two-pion exchange. These are the peripheral NN partial 
waves below the inelastic iViV7r-threshold. 

In ref. |l| the 27r-exchange contributions were worked out up to third order in small 
momenta including one-loop graphs with a vertex from the second order chiral ttN- 

(2) 

Lagrangian C w ^. This part of the Lagrangian involves several additional low-energy 
constants Cj, (i = 1,2,3,4), which have been recently determined from a fit to many 
low energy pion-nucleon data ||. Since some of these constants q are much larger than 
the natural scale 1/2M (M being the nucleon mass) they indeed produce the major 27r- 
exchange effect in NN-scattering. In particular, the isoscalar central potential is almost 
entirely given by the constant C3 proportional to the so-called nucleon axial polarizability 
As shown in ref. f2| it is the A(1232)-resonance that makes the dominant contribu- 
tion to the low energy constants 02,3,4 if one interprets their values in terms of resonance 
exchange. Therefore the 27r-exchange effects proportional to 03,4 as calculated in ref. 0] 
can be approximately identified as single A-excitation graphs for which one has ignored 
the energy dependence of the A-propagator by using just a contact ttttNN- vertex. How- 
ever, at the inelastic iViV7r-threshold the nucleon kinetic energy T^b becomes comparable 
to the AiV-mass splitting of 293 MeV, and a treatment of important A-dynamics via 
contact interactions may be too crude. This is also indicated by the results in ref. 
which show too large attraction in the F- waves above Ti a b = 180 MeV and in the D- waves 
even for smaller energies around 50 MeV. One can expect that explicit A (1232) degrees 
of freedom cure this problem (at least partly). 

In the past the A-excitation processes in NN-scattering have been considered in various 
approaches, either via dispersion theoretical methods |4]] or using iVA-transition potentials 
which neglect some (technically complicated) parts of the diagrams. At present no 
proper quantum field theoretical evaluation of these diagrams exists, which would avoid 
artificial "form factors" or cutoffs. As in ref. || such "form factors" are often introduced 
in order to enforce convergent loop integrals. In this work we evaluate the single and 
double A-excitation graphs in the static limit M — ► 00 (and also their first relativistic 
1/M-correction) using covariant perturbation theory and dimensional regularization. In 
one-loop order to which we are working here the divergences of the diagrams show up 
only as purely polynomial NN-amplitudes which do not contribute to the phase shifts 
with orbital angular momentum L > 2 and the mixing angles with J > 2. Therefore we 
can study the effects due to A-excitation in a completely parameterfree fashion. 

Next we consider, at two-loop order, the two-pion exchange diagrams involving the 
chiral 7T7r-interaction, so-called correlated 27r-exchange. In phenomenological approaches 
H^] such processes are often identified with the enhancement around masses of 550 MeV 
observed in the isoscalar central (tttt — > NN) spectral function proportional to |/o+| 2 - 
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Such a behavior is found in the dispersion theoretical analyses of the 7riV-scattering data 
and often called "a-meson". After transforming the respective NN-amplitudes into 
a local coordinate space potential, we find that these two-loop diagrams involving the 
chiral 7T7r-interaction lead to rather different effects. In particular the resulting isoscalar 
central potential turns out to be weakly repulsive, contrary to common believe. The 
repulsive nature of isoscalar central potential finds its explanation in the isospin-zero S- 
wave 7T7r-amplitude which becomes repulsive sufficiently far below the threshold. In any 
case, all effects due to correlated 27r-exchange (at two-loop) that we find here in (heavy 
baryon) chiral perturbation theory are in fact negligibly small. This means that the 
enhancement showing up in the spectral function \fo+\ 2 cannot simply be identified with 
the (perturbative two-loop) correlated 27r-exchange diagram. 

With decreasing orbital angular momentum L we find that the description of the 
empirical F-wave NN phase shifts above = 150 MeV requires a further (well-known) 
ingredient, namely the exchange of vector mesons, p and u. While tu-meson exchange 
(using g^ N /4:Tr = 12.9 for its coupling constant) provides the necessary overall repulsion, 
the p-meson with its large tensor-to-vector coupling ratio (k p = 6) leads to the correct 
splitting of the singlet and the three triplet F-waves. 

In D-waves we find deviations from the empirical NN phase shifts already at = 100 
MeV (partly even for lower T/ &)- F° r these low angular momentum partial waves the van 
der Waals behavior of the 27r-exchange NN-potential becomes problematic, since its r -6 - 
singularity practically extinguishes the natural centrifugal barrier effects coming from the 
wave function r 2 j2(pr) 2 (with j2{w) a spherical Bessel function). At such kinematics the 
short range NN-repulsion starts to become essential. Its dynamical origin lies of course 
outside the perturbative interaction of point-like baryons and mesons treated here. For 
some recent attempts on this problem in the framework of effective field theory see ref. 
MI 

We can summarize our work as follows: Within systematic perturbation theory based 
on chiral symmetry one finds an accurate description of the empirical NN phase shifts in 
the partial waves with L > 3 up to = 350 MeV and partly up to 80 MeV in D-waves 
with the following ingredients: 

1) Point-like one-pion exchange (without introducing an unmeasurable ttN form factor) 

2) Iterated one-pion exchange 

3) Irreducible two-pion exchange with only nucleons in intermediate states (plus their 
first relativistic 1/M-correction) 

4) Two-pion exchange with single and double A-isobar excitation (in the static limit) 

5) Vector meson (p and uS) exchange with standard values of the coupling constants 
but no ad hoc form factors 

These components comprise the dynamics of the peripheral nucleon-nucleon interac- 
tion. The emerging physical picture is of course not entirely new. However, our calculation 
is based on an effective chiral Lagrangian and we apply rigorous methods of perturbative 
quantum field theory (covariant Feynman graphs and dimensional regularization) . Thus 
there is no need to introduce extra cut offs or ad hoc form factors which sometimes obscure 
the real physics. 
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On the experimental side there are upcoming precision data from the Indiana Cooler 
Synchrotron Facility (IUCF), which will lead to an improved NN phase shift analysis in 
the energy range below the pion production threshold. We propose to use the chiral NN 
phase shifts with L > 3 presented here as input in a future phase shift analysis. 

II. BASIC FORMALISM 

In this section, we briefly review some basic formalism needed to describe elastic 
nucleon-nucleon scattering. In the center of mass frame the on-shell T-matrix for the 
process N(p) + N(—p) — > N(p') + N(-p') takes the following general form, 

Tnn = Vc + ti- f 2 W c + [V s + fx ■ f 2 W s ] 3\ ■ a 2 + [V T + fx ■ f 2 W T ] Bx-qo 2 -q 
+ [V so + n ■ t 2 W so ] + <?2) ■ (qxp) 

+ [Vq + fx ■ f 2 W Q ] o\ • (q x p) a 2 -(qxp). (1) 

The ten complex functions Vc, ■ ■ ■ , Wq depend on the center of mass momentum p = 
\p\ = \p' | and the momentum transfer q = \q \ with q = p' — p. The subscripts refer 
to the central, spin-spin, tensor, spin-orbit and quadratic spin-orbit components, each of 
which occurs in an isoscalar (V) and an isovector (W) version. In terms of the nucleon 



laboratory kinetic energy the center of mass momentum p is given as p = y T} a bM/2. 

In order to compute the NN phase shifts and mixing angles one needs the matrix 
elements of T N n in the LS J-basis, where L = J — S, J, J + S, S = 0, 1 and J = 0, 1, 2, . . . 
denote the orbital angular momentum, the total spin and the total angular momentum 
respectively. The explicit form of the partial wave projection formulas can be found in 
section 3 of ref . [Jxj] . Phase shifts and mixing angles are then given perturbatively as 

hsj = ^ Re {LSJ\T NN \LSJ), (2) 

ej=^Re{J-l,U\T NN \J + l,lJ), (3) 

with the center of mass nucleon energy E = \/M 2 + p 2 . The perturbative expressions in 
eqs.(2,3) apply only if the phase shifts 5lsj and mixing angles ej are sufficiently small in 
order not to have any substantial violations of unitarity. We calculate reliably only those 
phase shifts and mixing angles which are smaller than 10° in magnitude. For these cases 
perturbation theory is well justified. 



III. ONE LOOP DIAGRAMS WITH DELTA-EXCITATION 

The one-loop diagrams with single and double A-isobar excitation are shown in Fig.l. 
The first diagram of triangle shape is specific for a calculation based on a chiral effective 
7riV-Lagrangian, since it involves the Weinberg- Tomozawa isovector iViV7r7T-contact ver- 
tex. Since all Feynman rules for the nucleon vertices and propagator are well documented 
in the recent review ||12||, we give here only some details associated with the A-isobar. In 



the heavy mass limit |T3| where one considers the external momenta and the delta-nucleon 
mass splitting A = Ma — M = 293 MeV small compared to the nucleon mass, M = 939 
MeV, one gets for the delta-propagator and the A — > 7i a N transition vertex 
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k - A + iO+ 



^S-lT a 
2V2f n 



(4) 



Here k^ denotes the four-momentum of the propagating delta (with its energy k counted 
modulo the large nucleon mass M) and f is the four-momentum of the emitted pion with 
isospin a. The 2x4 spin and isospin transition matrices Si and T a satisfy the relations 



Si S 



(25, 



It; 



jk(Tk)/3 and T a T b T = (25 ab - ie abc T c )/3 Q- For the 7riVA coupling 
constant we have already inserted the large- N c value g n NA = 3g w N/\/2 together with 
the Goldberger-Treiman relation g^isi = QaM/ f n = 13.4 (with f n = 92.4 MeV this gives 
Qa = 1.32). For the A — > ttN decay width one finds in this case 



T(A -> ttN) 



3g\ E 



N 



M 



32vr/2 M + A 



(E 2 N - M 



2\3/2 



110.6 MeV, 



(5) 



a number which is in good agreement with the empirical decay width r(A — > nN) = 
(115 ± 5) MeV. Here E N = M + (A 2 - m 2 ) /2M A = 966 MeV denotes the center of mass 
energy of the decay nucleon. We note that eq.(5) for the A-decay width derives from a 



relativistic calculation using Rarita-Schwinger spinors 



Since the momentum of the 



decay pion is not so small, namely 227.3 MeV, one would loose important kinematical 



factors in a non-relativistic approximation and overestimate the A-decay width [H . 
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s ✓ 
\ ✓ 



Fig.l: One-loop 2n-exchange diagrams with single and double A(1232) -excitation 



Now we give the one-loop NN-amplitudes which result from evaluating the diagrams 
shown in Fig.l in the heavy mass limit, M — > oo. For the sake of simplicity we omit 
purely polynomial terms which do not contribute to the phase shifts with L > 2 and 
to mixing angles with J > 2 (see section 4.1 in ref. Jl[). The divergences of the loop 
diagrams are actually included in such (irrelevant) polynomial terms. We note that all 
one-loop vertex and self energy corrections with A-isobar excitation to the l7r-exchange 
diagram lead only to mass and coupling constant renormalization; they do not introduce 
a pion-nucleon "form factor" (see section 4.1 in [|IJ). We find the following analytical 
results for the three classes of diagrams shown in Fig.l: 



a) A-excitation in triangle graphs: 



W, 



c 



9 2 A 



192tt 2 /4 

b) Single A-excitation in box graphs: 



6S - w 2 )L(q) + 12A 2 S D(q)} . 



(6) 



c 



% ' -(2m 2 +g 2 ) 2 A(g) 



327T/4A 



(7) 



W, 



c 



192vr 2 /4 



|(12A 2 - 20m 2 - llq 2 )L(q) + 6S 2 D(q)} 
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V T = ~V S = ~ 2L(«) + (» 2 - 4A 2 )B(,)} , (9) 

"' T = 7 , ^124a" MW ' <10) 
c) Double A-excitation in box graphs: 

Vc = ~ 4A2L(g) + S ™ + (S + 8A2 ) D ^)]} ' (") 

Wc = ^Jl{^ - ^) L (q) + 3 s (?) + (§ A2 " Z)D(q)]} , (12) 

Vt = ~ V S = ^Jl{^{q) + (12A 2 - w 2 )D(q)} , (13) 

W T = -j 2 W s = I5 A^{2L(,) + (4A 2 + W 2 )D( g )} . (14) 

The following set of loop functions and abbreviations has been used to express these 
NN-amplitudes, 

L(q) = ^ In , w = y/Aml + q 2 , (15) 

A(q) = — arctan - , (16) 

w 2q 2m w ' v ; 



1 roo dfx y/^ 2 -4m2 
= A jL arCtan 2A ' (17) 



d/f _______ y/i 2 -4m^ 

2 m7r /i z + g 2 
2S 

S = 2m 2 +g 2 -2A 2 . (19) 



The isoscalar central Vc and isovector tensor amplitude Wt coming from the box 
graphs with single A-excitation in eqs.(7,10) show an interesting feature. The dependence 
on the delta-nucleon mass splitting A is a trivial factor A -1 . This means that, in the sum 
of the corresponding planar and crossed box graphs, the energy dependence of the delta- 
propagator has effectively disappeared. A 7T7riViV-contact vertex (— c 3 = 2c 4 = g\/2A) 
proportional to the inverse mass splitting A" 1 instead of the propagating A-isobars would 
give exactly the same result. Another way to understand this coincidence is to compute 
the sum of the energy denominators for all possible time-orderings which gives 

2oj 1 uj 2 (u; 1 + oj 2 ) + A(w 2 + 3luiLU 2 + uj 2 ) + A 2 (cji + u 2 ) 
Acjicj 2 (^i + u) 2 )(ui + A)(uo 2 + A) 

for the planar box graph and 

Ldj + UOiU) 2 + UJ 2 + A(u>i + LU 2 ) 



ujxuj^uji + uj 2 ){uji + A)(u 2 + A) 



(21) 



for the crossed box graph, with uj± and uj 2 denoting the on-shell energies of the two 
exchanged pions. The sum of the two expressions in eq.(20) and eq.(21) leads to the 
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surprisingly simple expression 2(Au;iCo>2)~ 1 m which the mass splitting A has factored 
out. In all other cases there is however a nontrivial dependence of the NN- amplitudes on 
the mass splitting A, in particular due to the occurrence of the loop function D(q) for 
which we have given its spectral representation in eq.(17). 

The NN-amplitudes in momentum space given above can be transformed into local 
coordinate space potentials (disregarding zero-range <5 3 (r )-terms) in the form of contin- 
uous superpositions of Yukawa functions (for details see section 6 in ref. |]). The mass 
spectra entering in this representation are given by imaginary parts of the NN-amplitudes 
analytically continued to time-like momentum transfer q = ifi — + . This requires the 
knowledge of the imaginary parts of the loop functions for (i > 2m 7n 



> 2 -K 

L{i(i) 



(i 

am = -L 



fi+ J(i 2 -Ami .vr 

In % — 

2m n 2 



(i + 2m 7T . 

In h 1 7r 

jj, — 2m n 



(22) 
(23) 



7T J (i 2 — Ami 

ImD(iii) = — arctan . (24) 

v p; 2(i A 2A v ! 

Here we have also given the real parts of L(i(i) and A(i(i) which will be needed later in 
the discussion of the correlated two-pion exchange. As in our previous work |1| we will 
denote r-space potentials by a tilde, i.e. Vb(r) is the isoscalar central potential, etc. 

The attractive isoscalar central potential Vc(r) generated by the A-excitation diagrams 
is shown in Fig. 2. The main contribution comes from the four graphs with single A- 
excitation and it has the following simple analytical form 

o 4 — 2x 

^" a,(r) = "mHa V (6 + 121 + 1012 + 4 * S + xi) 

= (2[W« 1 "(r)] 2 + m"\r)f} , (25) 

with the abbreviation x = m n r. We have expressed this potential in the second line in 
terms of squares of the l7r-exchange isovector spin-spin and tensor potentials Wjy (r). 
For the isovector spin-spin and tensor potentials coming from the single A-excitation 
graphs one finds, similarly, 

4 —2x 

= ^{[W^\r)] 2 - [Wt\r)] 2 } , (26) 

4 — 2x 

^ NA) (0 = -I9^%A V (1 + X)(3 + 3X + ^ 

= ^W^\r){wt\r) ~ W$*\r)} . (27) 

The relative weights between the l7r-exchange spin-spin and tensor potential in 
eqs. (25,26,27) can be easily understood from relations between the tensor operator 
5*12 (r) = 3<t ■ f c?2 • f — c?i • &2 and the spin-spin operator a\- fii- These relations are: 
[S 12 (f)} 2 = Q + 2a 1 -a 2 - 2S 12 (f), {S 12 (f), a l ■ a 2 } = 2S 12 (f) and {B x ■ a 2 f = 3 - 2a x ■ a 2 . 
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FIGURES 




Note that the relevant energy difference A/12 = 24.4 MeV in eq.(25) is rather small. 
It is quite remarkable that these simple relations between Vq NA, \t), (r) and the 

l7r-exchange potentials (eq. (25,26,27)) are exact for the isoscalar central and the isovector 
spin-spin/tensor potentials generated by the direct and crossed box graph with single A- 
excitation (in the static limit M — > oo). Heuristic considerations already suggest such a 
form of the isoscalar central 27r-exchange potential. 

The double A-excitation diagrams amount to about 30% of the total isoscalar cen- 
tral potential as shown by the dashed-dotted line in Fig.2. The fictitious 'V'-exchange 
potential % a \r) = -{g 2 j 1 ^r)e~ M ° r with gl/Ax = 7.1 and M a = 550 MeV |§ is also 
shown in Fig.2 by the dotted line. For distances r > 2 fm one observes an almost perfect 
agreement between the phenomenological "a" -exchange potential and the total isoscalar 
central potential generated by 27r-exchange with A-excitation. For shorter distances r < 2 
fm the latter one increases more strongly in magnitude due to its inherent r _6 -singularity 
(van der Waals behavior). Among all potentials generated by A-excitation the isoscalar 
central Vc(r) is by far the largest one. The attractive isovector central potential Wc{t) 
is approximately a factor ten smaller, with values of —32.7 MeV, —2.0 MeV and —0.25 
MeV at distances r = 1 fm, 1.5 fm and 2.0 fm. The repulsive isoscalar and isovector 
spin-spin potentials Vs(r), Wsir) and the attractive isoscalar and isovector tensor poten- 
tials Vt(t), Wr(r) are even smaller, with typical values of ±13 MeV at r = 1 fm. The 
asymptotic behavior for large r is e - 2m ^ r r - 5 / 2 f or the isovector central potential Wc(r) 
and e - 2m ^ r r - 7 / 2 f or the isoscalar spin-spin/tensor potential Vs,T(f)- In the other cases it 
can be read off from eqs. (25,26,27). We note that throughout the large-r asymptotics is 
determined alone by the box graphs with single A-excitation b). 

We have furthermore evaluated the first relativistic correction, proportional to 1/M, 
to the A-excitation graphs starting from the Rarita-Schwinger form of the spin-3/2 prop- 
agator and 7riVA-vertex. Explicit formulas for the respective NN-amplitudes can be found 
in the appendix. Unfortunately these corrections are not small in the range 1 fm < r < 2 
fm or for Ti a b > 100 MeV. The reason is probably found in combinatoric factors (up to six 
vertices and propagators are expanded in powers of 1/M) and in the fact that the ratio 
A/M ~ 0.3 is not so small. We expect that part of these relativistic 1/M-corrections will 
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be cancelled by higher orders in the 1/M-expansion. On the other hand it is not clear 
how realistic a description of the A-resonance via a local Rarita-Schwinger spinor-field 
really is, even though it is the only available Lorentz-covariant formalism for spin-3/2 
particles |nj. Off its mass-shell the Rarita-Schwinger spinor-field exhibits potentially un- 
physical spin-1/2 degrees of freedom. In view of these problems we prefer not to include 
the relativistic 1/M-correction to the A-excitation graphs (Fig.l) in our calculation. 

Let us give another argument in support of the idea that the A-excitation in the 
static limit (M — ► oo) already represents the low energy titt — ► NN dynamics in the 
scalar isoscalar channel quite well. The shift of the nucleon scalar form factor crjv(i) 
from momentum transfer t — to t — 2m 2 has been evaluated in ref. |L7| via dispersion 
relations and empirical tcN- and 7T7r-data. The result is aN(2m 2 ) — o"tv(0) = 15.2 ± 0.4 
MeV. In the one-loop approximation of heavy baryon chiral perturbation theory one finds 
for this scalar- isoscalar quantity |fL8| , 



^nlm 2 ( i A + \/~K~ 2 

a N (2ml) - a N (0) = gg^fm. - 4^^ ln V_ 



ml 



64tt 2 /2 

+ (tt - 4)A + 8A 3 D(iV2m n ) \ = (8.0 + 6.9) MeV (2* 



where the first term ~ irm n comes from the pion loop-diagram with a nucleon intermediate 
state, and the remaining ones come from an analogous diagram with an intermediate 
A(1232)-isobar (8A 2 D(iV2m n ) = 7.02). The sum 14.9 MeV of both terms agrees well 
with the dispersion-theoretical value 15.2 ±0.4 MeV. While the detailed spectral function 
Im<7Ar(t) given by these two one-loop diagrams is not in perfect agreement [Tj| with the 



empirical one derived via dispersion theory Hl7 |, the relevant integral over lvsia^it) /t{t — 
2m 2 ) is well reproduced. 

We expect a similar mechanism to be at work for the 27r-exchange isoscalar central NN- 
potential which can anyhow only be tested at large and intermediate distances, r > 1.5 
fm. The mass spectrum ImVc as given by the A-excitation diagrams in the static limit 
does certainly not have all structures as the one derived from dispersion theory j|,|9]], in 
particular there is no enhancement around 550 MeV (the broad "a" -meson). By analogy 
with the previous discussion of the nucleon scalar form factor we expect that the long and 
intermediate range components of the isoscalar central potential Vc(r) are well represented 
by the A-excitation graphs in the static limit. For the long and intermediate range 
isoscalar central potential only the low energy part and some global features, but not the 
details of its mass spectrum, play a role. The good agreement between the "a" -potential 
and the A-excitation 27r-exchange potential for r > 2 fm in Fig. 2 demonstrates this fact 
very clearly. 



IV. CORRELATED TWO-PION EXCHANGE 

As previously mentioned the mass spectrum in the isoscalar central channel ImVc ~ 
|/o+| 2 , with the tttt — > NN S-wave amplitude /o+, shows an enhancement around 550 
MeV (the so-called "a" -meson) 0,0. This enhancement is often interpreted as originating 
from correlated two-pion exchange (the first diagram shown in Fig. 3) in which the two 
exchanged pions interact while propagating between the two nucleons. This motivated us 
to evaluate this specific two-loop diagram in heavy baryon chiral perturbation theory. 
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Fig. 3: Correlated two-pion exchange diagrams and some related two-loop graphs. The 

dot represents the tttt -interaction 



The correlated 27r-exchange diagram involves the off-shell 7T7r-interaction which is, 
however, not uniquely defined by the effective chiral Lagrangian. It depends on the choice 
of the interpolating pion field or, to be explicit, on the parametrization of the SU(2)-matrix 
U(tt ) which is the basic variable entering the non-linear chiral 7T7r-Lagrangian, 



1 + TJ 



71 



-7T ' 



% a 

7*' 
j - 



T ■ 7T 



a — 



f 1 

J 7T 



■-TT*+... 



(29) 



The arbitrariness of the coefficient a in the 3rd and 4th power of the pion field just reflects 
this ambiguity For the pure, isolated 7T7r-interaction this does not cause a problem since 



the off-shell amplitude A Tl 



-2a (4ml— Qi~ Q2~ qD)/ f% * s n °t an observable 



and the on-shell amplitude can be easily shown to be independent of a. However, the 
isoscalar central NN-amplitude Vc generated by the first diagram in Fig. 3 is a-dependent. 
It is therefore not meaningful to consider this diagram alone, instead one has to find the 
complete subclass of diagrams for which the unphysical a-dependence drops out. This 
subclass is obtained by adding those diagrams which result from shifting the 47r-vertex 
to a nucleon line (see a typical representative in Fig.3). This procedure resembles the 
construction of gauge invariant subclasses. The four additional diagrams now indeed 
cancel an unwanted term proportional to 10a — 1. We note that for isovector amplitudes 
there is no a-dependence in the first diagram of Fig.3. Altogether we find the following 
contribution to the NN T-matrix from the factorizable two-loop diagrams in Fig.3: 



V c 
W c 



1024tt 2 ^ 
1 



ml + 2q 2 ) [m ff + {2ml + Q 2 ) A (q) 
Aml(l + 2g 2 A ) + q 2 (l + 5g 2 A )}L(q) 

2 2 



(30) 



184327r 4 /^ 

Aml(l + 2g\) + q 2 (l + 5g 2 A ) In ^ - |(5 + 13^)} 



W 7 



--W s 
q 2 



20487r 2 /^ 

*> TV 



m w + w 2 A(q) 



(31) 
(32) 



Note that Vc and Wt are finite in dimensional regularization. The divergent terms in 
Wc have been omitted according to the usual minimal subtraction prescription which 
introduces a logarithmic dependence on the renormalization scale A. The expressions in 
eqs. (30,31,32) for correlated 27r-exchange in heavy baryon chiral perturbation theory have 
a simple interpretation in terms of the one-loop contributions to the nucleon scalar and 
isovector electromagnetic form factors (19| : 
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V c 
W c 
W T 



1 



8/ 2 
1 

~~2 



9 



W s 



G\)(—q 2 )ioop 
1 



32M 2 / 2 



9 )loop 



(33) 
(34) 
(35) 



With to( S 7T7r 



(25,,^ — m 2 )/(327r/ 2 ) the on-mass-shell isospin-zero S-wave 7T7r-amplitude 
at leading order. Actually, we have considered for Wc a larger class of factorizable two- 
loop graphs as shown in Fig. 3, such that no additive constant to the electric from factor 
<? 2 )ioo P appears under the square in eq.(34). Note that both central NN-amplitudes in 
eqs. (33,34) are negative which indicates a repulsive NN-interaction. The negative sign of 
Vc comes from the isospin-zero S-wave 7T7r-amplitude evaluated at negative s nn = —q 2 < 0. 
As a matter of fact chiral soft pion theorems require the 7T7r-interaction to be of derivative- 
nature which unavoidably leads to strong energy dependence of the isospin-zero S-wave 
amplitude. An immediate consequence is that the attractive isospin-zero S-wave tttt- 
interaction above threshold (s W7r > 4m 2 ) switches over to repulsion sufficiently far below 
the threshold (s W7r < 0). In the NN-scattering process the exchanged pion pair belongs 
exclusively to the latter kinematical region < 0. This is the basic reason why we find 
here a (weakly) repulsive isoscalar central potential from correlated 27r-exchange. 



r[fm] 


1.0 


1.5 


2.0 


r - 


-> 


r — ► oo 


V c 


0.731 


0.090 


0.016 


r 


-4 


g 2 771 77 T ty, 2 Jj-^ 


W c 


3.811 


0.571 


0.092 




lnr 


e -2m 7r r r -5/2 


W s 


-0.170 


-0.022 


-0.004 


r" 


-4 


g 2777-7T V j, 3 


W T 


0.193 


0.022 


0.004 


r" 


-4 


g 2 777 Tj- T 3 



Tab.l: Coordinate space potentials from correlated In-exchange in units of MeV. 



In Tab.l we display some values for the coordinate space potentials as they derive 
from our calculation of the correlated 27r-exchange using the corresponding imaginary 
parts (Im [A(iyu)] 2 = 2KeA(i/j,)lmA(i/j,) etc.) to evaluate their spectral representations 
|Ij]. In order to get an estimate of the magnitude of Wc(r) we set A = M w = 782 MeV. 
Indeed both central potentials and the tensor potential are repulsive but they are also 
negligibly small. We have confirmed this furthermore by evaluating the contribution of 
the NN-amplitudes in eqs. (30,31,32) to the phase shifts and mixing angles. In the case of 
the correlated 27r-exchange isovector spin-spin potential one can even give an expression 
in closed form, 

4 — 2x 

Ws{r) = 614 ^ 3/6 e r7 {(15 + 30a; + 24x 2 + 8x 3 )( lE + ln4x) 

+ (15 - 30x + 24x 2 - 8x 3 )e 4x E 1 (Ax) - 4x(15 + I5x + 8x 2 + 2x 3 )} , (36) 

with x = m n r, 7^ = 0.5772... the Euler-Mascheroni number and Ei(Ax) = f££ d££ -1 e~^ 
the exponential integral function. The isovector spin-spin potential Ws(r) in eq.(36) 
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has only a r _4 -singularity near the origin and it vanishes identically in the chiral limit, 
= 0. The isoscalar central potential Vc{r) and the isovector tensor potential Wt(t) 
from correlated 27r-exchange have actually the same features. In contrast to this the 
isovector central potential Wc{r) from correlated 27r-exchange does not vanish in the 
chiral limit m n = and it has a much stronger r~ 7 In r-singularity near the origin. In 
Tab.l the asymptotic behavior for large r is also given. 

In ref. || further correlated 27r-exchange diagrams with single and double A-excitation 
were considered. We restrict ourselves here to a rough estimate of their effects. Obvi- 
ously the contribution of the diagrams with A-excitation is overestimated in the limit of 
zero mass splitting, A = 0. In that case only different spin/isospin factors show up in 
comparison with the diagrams involving nucleon intermediate states. The isoscalar cen- 
tral amplitude Vq in eq.(30) would simply be multiplied by a factor 9. Comparing with 
the scalar form factor difference <7jv(2mJ;) — <7jv(0) where the A-isobar in the zero mass 
splitting limit triples the pure nucleon-term, one expects an amplification factor around 
4 to be a more realistic estimate of A-effects for finite mass splitting. In the case of the 
isovector tensor amplitude Wt the degenerate A-isobars would lead to a factor 9/4 in 
eq.(32), and for the isovector central amplitude Wo the ^-dependent terms in eq.(31) 
would be cancelled by degenerate A-isobars. Even including such amplification factors 
due to intermediate A-isobars the correlated 27r-exchange is a negligibly small correction. 
This completes the discussion of the correlated 27r-exchange in heavy baryon chiral per- 
turbation theory. As a main result we find that this diagram cannot be identified with 
the enhancement showing up in the empirical isoscalar central spectral function |/b+| 2 - 
Analogous features have already been observed in the one-loop calculation of the ttN 
scattering amplitude in ref. ||. There the (class of) diagrams with pion self-interaction 
reduces the spin- and isospin averaged P-wave scattering volume Pi as well as the nucleon 
axial polarizability a a- A reduction of these quantities is just the opposite of what one 
naively expects from scalar-isoscalar 7T7r-correlations. 

V. VECTOR MESON EXCHANGE 

It is clear that the coupling of vector mesons (p, iS) to nucleons has no direct relation to 
chiral symmetry. The vector meson mass scale actually marks the kinematical endpoint of 
chiral dynamics with weakly interacting pions. However, vector mesons are an important 
ingredient for the understanding of nucleon structure, as witnessed e.g. by the dispersion- 
theoretical analysis of the nucleon electromagnetic form factors P,^0| and the successful 
one-boson-exchange models of the NN-interaction || . In our study of the peripheral NN 
phases we find that the chiral 27r-exchange considered so far is insufficient to describe the 
empirical F-wave phase shifts above TJ a & = 150 MeV (see Fig.4 in our previous work [0). 
The source of the discrepancy is a much too strong attraction between the nucleons for 
distances r < 1.5 fm. Naturally one expects that vector meson exchange produces the 
repulsive interaction at intermediate distances which will be able to cure this problem. 

The coupling of the p- and to- meson to the nucleon is characterized by vector cou- 
pling constants g P N,uN an d tensor-to- vector coupling ratios k Pj£J . We will apply symmetry 
relations to minimize the number of free vector meson parameters. The KSFR relation 
M p = = V^gpf-K = 782 MeV leads to g p = 6.0 for the universal p-coupling constant. 
Via isospin considerations one obtains g p ^ = g p /2 = 3.0 and SU(3) symmetry would pre- 
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diet g 



SU(3) 



3g p /2. However, all existing determinations of g^N point towards a larger 

SU(3) 



jN 

value. We use here g^N = y/^g^^' = 12.7 as obtained from a global fit of the NN F 



waves. This number is in the range of values found from forward NN-dispersion relations 
Ull (duiN = 10.1 ± 0.9), the Bonn NN-potential [§] (g^N = 11.5) and the recent disper- 
sion analysis of the nucleon electromagnetic form factors [PJ (g wN = 20.9 ± 0.3). For 
the tensor-to- vector coupling ratio we use the values of ref. ||, k p = 6 and k u = 0, 
which are also confirmed by the nucleon electromagnetic form factor analysis of [[2(J 
(k p = 6.1 ± 0.2, k w = —0.16 ± 0.02). The precise value of is irrelevant for the NN 
F-wave phase shifts as long as it is sufficiently small, \k u \ < 0.2. The p-meson exchange 
between two nucleons leads to the following isovector NN-amplitudes, 



a 2 



~ AM 2 (M 2 p +q< 



o V F + M MJ 



+Up 2 -q 2 )( ^ 1 



2^E + M M>\ 

2 



W T = -\w s 



AM 2 (M 2 + q 2 ) ' 



(37) 
(38) 



Wso = 



n 2 

y P N 



~4M 2 {M 2 + q 2 
E 



2E- q - 



2^E + M M 



+ -)}( 



i 



E + M M> 



1 + kJ— - 



Wn = 



W AM(E + M) 



K p q* 



2M(E + M) 

9p N \ (V-g 2 )^ / 1 Kn \ - 

AM 2 {M 2 + q 2 ) \4:M 2 (E + M) 2 ^E + M 



(39) 



• 2 V 



M> 



2K f 



' M(E + M) 



(E q 2 
Kp W AM(E + M) 



(40) 



where we have kept the fully relativistic expressions (E = ^M 2 + p 2 ). In a non-relativistic 
truncation at order M~ 2 one would otherwise loose important contributions, in particular 
the quadratic spin-orbit term Wq proportional to the large coefficient 2n 2 p + 2k p + 1/4. 
Interestingly, the fully relativistic tensor and spin-spin terms Wt,s m eq.(38) agree exactly 
with the lowest order non-relativistic approximation. For a pseudoscalar meson exchange 
this is also the case [l]. The isoscalar NN-amplitudes due to u-exchange are obtained 
in complete analogy by the replacement (g P N,K p ) — > (g W N,Ku) in eqs. (37-40). For u- 
exchange (with ~ 0) the truncation at order M~ 2 is sufficiently accurate. At energies 
where p- and tu-exchange are important it is not meaningful to approximate them by local 
(polynomial) contact terms, since the ratio q/M w is not small. 

For the sake of completeness one should also add the exchange of ^-mesons with mass 
m v = 547.45 MeV. In the absence of a reliable empirical determination of the ?7iV-coupling 
constant g V N we will use the SU(3)-value g V N = (3F — D)M/ VSf^ = 4.4 together with the 
approximate values of the octet axial vector coupling constants D = 3/4 and F — 1/2. For 
comparison the Bonn OBE-model (without 27r-exchange) uses g V N = 6.8. The ^-exchange 
leads to an isoscalar tensor amplitude of the form, 



V r 



64/2( m 2 +? 2) 



(41) 
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The actual calculation of the phase shifts shows that //-exchange with the coupling 
strength given by SU(3) is almost negligible in all partial waves with L > 2. We have fur- 
thermore evaluated the K K -exchange (bubble and triangle) diagrams in heavy baryon chi- 
ral perturbation theory. As expected these processes lead to very small repulsive isoscalar 
and isovector central NN-potentials, V c {lfm) = 0.99 MeV, W c {lfm) = 0.26 MeV. 



VI. RESULTS FOR PHASE SHIFTS AND MIXING ANGLES 



In this section we present and discuss our results for the NN phase-shifts with L > 3 
and mixing angles with J > 3 up to nucleon laboratory kinetic energies of Ti ab = 350 MeV. 
For the D-wave phase shifts and €2 we show results only up to T/ a b = 120 MeV. First, we 
state all ingredients which go into the calculation. We include all terms derived in the 
recent work |IJ, the point-like l7r-exchange, the iterated l7r-exchange and the irreducible 
27r-exchange with the low energy constants ci^ 4 set equal to zero. The terms proportional 
to c 3i4 are now substituted by explicit A(1232)-dynamics, and c\ gave anyhow only a 
marginal contribution to the isoscalar central amplitude Vc- The new ingredients are 
the 27r-exchange with A-excitation in the static limit (section 3) and the vector meson 
exchange (section 5) with coupling constants g P N = 3, k p = 6, g^N = 12.7, = 0. We 
use throughout the value = 138 MeV for the (average) pion mass. 



A. D- WAVES 



The D-wave phase shifts and mixing angle 62 are shown in Fig. 4 up to TJ^ = 120 
MeV. The dashed curve corresponds to the one-pion exchange approximation and the full 
curve includes in addition two-pion exchange and vector meson exchange. The dotted 
curve represents the recent empirical energy dependent NN phase shift analysis of ref. 
p2|| (VPI). The triangles and squares give NN phase shifts and mixing angles derived 
from single energy analyses of ref. |23| and ref. , respectively. The circles represent the 
results of the multi-energy partial wave analysis of ref. |24|] (Tab.IV,V). In all cases the 
two-pion and vector meson exchange corrections go into the right direction, but deviations 
show up already above Ti ab = 30 MeV in the 1 D2 partial wave and above T\ ab = 50 
MeV for 3 D\ and 3 D^. The 3 D2 phase shift and the mixing angle €2 are in agreement 
with the data up to Ti ab = 100 MeV. Similar results were found recently in ref. 
using a somewhat different approach to the 27r-exchange. Compared to our previous 
calculation [[[]], there is no improvement in the D-wave phase shifts and 62 due to adding 
(perturbative) p- and cj-exchange. The 27r-exchange with its r -6 singular behavior still 
provides too large attraction at distances r < 1 fm. Obviously, the D-wave phase shifts 
and 62 above Ti ab = 100 MeV are already sensitive to the short range NN-repulsion beyond 
cu-exchange. It appears, that the D-waves (above T[ ab = 100 MeV), as well as the S- and 
P-waves, require non-perturbative methods and phenomenological parametrizations of the 
short range NN-interaction. This is, of course, well known from earlier investigations. 



B. F- WAVES 



The F-wave phase shifts and the mixing angle e 3 are shown in Fig.5. We present here 
results up to Ti ab = 350 MeV, i.e. 70 MeV above the iViW-threshold where inelasticities 
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are still negligible. The phase shifts in the 1 F 3 , 3 F 2 , 3 F 3 partial waves are in very good 
agreement with the empirical data up to Ti a b = 350 MeV, whereas deviations show up 
in the 3 _F 4 partial wave above T} a b = 220 MeV. The good agreement in the former case 
results from the inclusion of p- and u-exchange. The cu-exchange compensates the too 
strong attraction from to 27r-exchange. The p-exchange with its large tensor, spin-orbit 
and quadratic spin-orbit amplitudes (k p = 6) leads to the correct splitting of the singlet 
and the three triplet F-wave phase shifts. In particular, the correct downward bending of 
the 3 i<2 phase shift is a vector meson exchange effect (see Fig. 4 in QX|] where the opposite 
behavior was found from 27r-exchange alone). The deviation in the 3 _F 4 partial wave above 
Tiab = 220 MeV suggests that further short range effects are at work in this particular 
channel. The mixing angle e 3 is in perfect agreement with the data for all energies up 
to Ti a b = 350 MeV. Irreducible 27r-exchange and iterated l7r-exchange contribute to this 
quantity with roughly equal strength. There is also a small but non-negligible contribution 
from vector meson exchange to e 3 above Ti a b = 200 MeV. 

C. G- WAVES 

The G-wave phase shifts and the mixing angle e 4 are shown in Fig. 6. The predictions 
are in good agreement with the data for all four partial wave phase shifts, and with the 
mixing angle up to = 350 MeV. The vector mesons p and uo almost do not operate 
anymore in these high angular momentum states, L = 4. At Ti a b = 300 MeV they produce 
phase shift contributions of 0.3° and smaller. One has now reached the chiral window in 
which In- and (chiral) 27r-exchange describe the NN-interaction completely and reliably. 
Note that the differences between l7r-exchange and data are still sizeable in the 1 G 4 and 
3 G5 partial waves. The chiral 27r-exchange closes this gap between the data and the 
l7r-exchange approximation. In the 1 G 4 partial wave the correction comes mainly from 
irreducible 27r-exchange, whereas in the 3 G§ partial wave (with total isospin 1 = 0) both 
irreducible 27r-exchange and iterated l7r-exchange contribute with roughly equal strength. 
Compared to our previous calculation the (very small) 3 G 5 phase shift does not bend 
over to positive values any more. For the 3 G3 and 3 G 4 phase shifts and e 4 the 27r-exchange 
corrections are relatively small. Nevertheless, these small effects improve the agreement 
between data and chiral predictions. 

D. H- WAVES 

The H-wave phase shifts and the mixing angle e§ are shown in Fig. 7. The 1 H$ phase 
shift and the mixing angle €5 have converged to the l7r-exchange approximation. In the 
3 if 4 and 3 if 5 partial wave one finds small corrections to l7r-exchange which nevertheless 
improve the agreement with the empirical data. Note however that the l7r-exchange 
approximation considerably underestimates the empirical 3 Hq phase shifts. Again this 
gap is closed by the chiral 27r-exchange. It is rather remarkable that the 27r-exchange 
effects are still important at such a large orbital angular momentum, L = 5. 
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E. I- WAVES 



The I-wave phase shifts and the mixing angle €q are shown in Fig.8. Again, the 3 Iq 
phase shift and the mixing angle have converged to the l7r-exchange. In the 
and 3 1? partial waves we predict small differences to the empirical phases of ref. ||22|| . The 
calculation presented here is indeed most reliable in the high angular momentum partial 
waves. 



F. INTERACTION DENSITIES 



In order to learn about the relevant length scales at which the peripheral NN inter- 
action actually takes place it is most instructive to study the local interaction densities. 
These are expressions of the form r 2 ji{pr) 2 ULSj{ r ) , with ]l{pt) a spherical Bessel func- 
tion and Ulsj( t ) the coordinate space potential in a given state \LSJ). Unfortunately, a 
fully equivalent representation of the phase shifts Slsj and the mixing angles ej in terms 
of local interaction densities is not possible because of some inherent non-localities of 
the iterated l7r-exchange (section 4.3 in ref. [|TJ) and the quadratic spin-orbit interaction. 
Nevertheless, the generic features are already displayed by the dominant contributions to 
the NN T-matrix, namely l7r-exchange and the isoscalar central components of 2ir- and 
cu-exchange. Two examples of such interaction densities are shown in Fig. 9 for the 3 _F 4 
partial wave at T[ ab = 200 MeV and the 1 G 4 partial wave at T [ab = 250 MeV. One can 
clearly see the peaking of the interaction densities around distances of about r = 1.5 fm 
and r = 2.5 fm, respectively. It results on one side from the centrifugal barrier effect 
given by the wave function r 2 ji{pr) 2 and on the other side from the exponential decay of 
the potential ULSj{ r )- Varying the orbital angular momentum L and the center of mass 
momentum p, the peak of the interaction density moves and appears approximately at a 
distance r « L/p. This length scale corresponds just to the classical impact parameter. 



VII. CONCLUDING REMARKS 



The chiral perturbation theory calculation of the NN phase shifts and mixing angles, 
presented here, is most reliable in the high angular momentum partial waves which probe 
the long and medium range parts of the NN force. This is the region where In- and 
27r-exchange explains the NN-interaction in a model independent way. We have demon- 
strated that the S- and P-wave chiral dynamics of the pion-nucleon system determines 
the peripheral NN phase shifts almost completely, with no arbitrary cutoffs or "form 
factors". The new aspect emphasized in this work is that the chiral pion-baryon effec- 
tive Lagrangian provides a well-defined systematic framework to deal with the peripheral 
nucleon-nucleon interaction. In particular there is no need to introduce the scalar-isoscalar 
"cr" -meson. The intermediate range isoscalar central attraction is explicitly produced by 
van der Waals-type 27r-exchange including intermediate A-isobar excitations. Effects 
from 7T7r-rescattering turn out to be negligibly small, in accordance with the suppression 
of higher loops in chiral perturbation theory. 
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Fig. 4-' D-wave NN phase shifts and mixing angle €2 versus the nucleon laboratory 
kinetic energy T^. The dashed curves correspond to the In -exchange approximation and 
the full curves include chiral 2tt- and vector meson exchange as well. The dotted curves 
represent the empirical energy dependent NN phase shift analysis of ref. j^jj . 
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Fig. 5: F-wave NN phase shifts and mixing angle €3 versus the nucleon laboratory 
kinetic energy Ti a b. For notations see Fig. 4. 
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Fig. 7: H-wave NN phase shifts and mixing angle 65 versus the nucleon laboratory 
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Fig. 8: I-wave NN phase shifts and mixing angle e§ versus the nucleon laboratory kinetic 
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Fig. 9: Examples of NN interaction densities in coordinate space. Dashed lines show 
the lir-exchange contributions. Full and dashed-dotted curves give the isoscalar central 
components of 2tt- and u- exchange, respectively. 
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APPENDIX: RELATIVISTIC 1/M-CORRECTION TO A-EXCITATION 

Here we will collect explicit formulas for the first relativistic correction proportional to 
1/M arising from the single and double A-excitation graphs in Fig.l. A convenient way 
to obtain these corrections is to use relativistic propagators and vertices and to perform 
the 1/M-expansion inside the loop integral. For the A(1232) with spin-3/2 this requires 
the Rarita-Schwinger formalism which gives for the A-propagator 

7 • k + Ma f 2k fl k u k^ 

3&«, - 7m7, - jTf2 + 77 ( 42 ) 



3P-M| + zO+l^ ;/ "'" Ml M L 
and for the A — > -K a N transition vertex 

3^7 



lV2f, y > 



(43) 



Z 1 is an off-shell parameter lying within in the empirically determined band —0.8 < Z' < 
0.3 |16|]. Again, we omit here additive linear polynomials cq 2 + c' in the central amplitudes 
and additive constants in the tensor and spin-orbit amplitudes which contain the divergent 
pieces of the diagrams. We will give separately the contributions coming from the three 
classes of A-excitation diagrams. We present first the results for Z' = and then display 
additional Z'-dependent terms as far as they show up at order 1/M. 

a) 1/M-correction to triangle graphs with A-excitation: 

W C = 192 f 2 M/4 { (4m - + 7q2)L{q) + 3Sg(g) + 3(4AV ~ ^ D ^} > ( 44 ) 
W so = 2W T = -j 2 W s = u ^ MfA { - 2L(q) + (w 2 - 4A 2 )D(g)} . (45) 
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b) 1/M-correction to box graphs with single A excitation: 

Vc = ^iW 2 ^ + 2q2K + ^ + 4< " 8mlw ~ 2]L{q) 

+ (2m 2 + q 2 ) 3 A-\A(q) + E [(2ml + <f ? ~ 4g 2 A 2 ]D(g)} , (46) 

= 3847r 4 /4A {3(2m 2 + 9 2 ) 3 A-VA( 9 ) + 2A 2 (8m 2 - q 2 )L{q) 

+3S 2 #(g) + 3S [(2m 2 + g 2 ) 2 + 4A 2 (g 2 - A 2 )]D(q} , (47) 

+E [2L(q) + (4A 2 + w')D(q)\} , (48) 

+2(4A 2 + E)L(g) + £ (u> 2 - 4A 2 ) J D(g)} , (49) 

+2S [2L(q) + (4A 2 + w 2 )£>(g)]} , (50) 

-2E L(g) + (4A 2 + S)(4A 2 - w 2 )D(q)} . (51) 
Additional contributions for Z' ^ 0: 



^ = -^Mf^ Z ' + 1)i2ml + q2fA{q) ' (52) 

Wt = -j 2 W s = ^L-QZ' + l) W 2 A( g ) . (53) 

c) 1/M-correction to box graphs with double A-excitation: 

+3[A- 2 S 3 + 8(g 4 + 6g 2 m 2 + 2g 2 A 2 + 8m 4 - 8m 2 A 2 )]D(g) 

+3[A- 2 (2m 2 + q 2 ) 2 + 4(2m 2 + 3g 2 - 3A 2 )]#(g)} , (54) 

Wc = 1536^M/^ 8(12A2 " 8ml + q2)L{q) ~ K{q) 

+3[16(2A 4 - 2ml + 3<? 2 A 2 - 3g 2 m 2 - q 4 ) - A~ 2 j: 3 }D(q) 

+3[4(6m 2 + 5g 2 - 5A 2 ) - A~ 2 (2m 2 + q 2 ) 2 ]H(q)} , (55) 

-E[2L(q) + (12A 2 + w 2 )D(q)}} , (56) 
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+8A 2 H(q) + (8A 2 + S)(w 2 - 4A 2 )D(q)} . 



(57) 



Additional contributions for Z' 



V c 




i 

— {3S [2A 2 — E — Z'(2A 2 + 3S)]D(g) 



+2[7m 2 + 4g 2 - 6A 2 + Z'(13m 2 + 7g 2 - 6A 2 )]L(g)} , 



(58) 



W T 




(59) 



The abbreviation K(q) stands for 



K(q) 



48 S 3 



L(q) 



w 2 - 4A 2 



=L'(2y/A*-ml)-L(2y/A*-ml) , (60) 



( w 2 _ 4A 2)2 



4^A 2 - m: 



TT 



and we refer to eqs. (15-19) for the definition of the functions w, S, L(q), A(q), D(q). Note 
that to order 1/M only the isoscalar central Vc and the isovector spin-spin/tensor Ws,t 
NN- amplitudes depend on the uncertain off-shell parameter Z' . Some terms given above 
have a prefactor g\j \M f^A 2 ) of dimension mass -7 . In coordinate space these terms will 
lead to a potential with a rather problematic r _8 -singularity near the origin. 

In particular, we have calculated here for the first time the spin-orbit interac- 
tion (V S o, Wso) generated by 27r-exchange with A-excitation. These spin-orbit NN- 
amplitudes are truly relativistic effects (absent in the static limit and independent of 
Z') and they may be of interest in future studies of the NN spin-orbit interaction. 
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